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[i] Aerosol-cloud interaction studies to date consider aerosol with a substantial fraction of 
soluble material as the sole source of cloud condensation nuclei (CCN). Emerging evidence 
suggests that mineral dust can act as good CCN through water adsorption onto the 
surface of particles. This study provides a first assessment of the contribution of insoluble 
dust to global CCN and cloud droplet number concentration (CDNC). Simulations are 
carried out with the NASA Global Modeling Initiative chemical transport model with an 
online aerosol simulation, considering emissions from fossil fuel, biomass burning, marine, 
and dust sources. CDNC is calculated online and explicitly considers the competition of 
soluble and insoluble CCN for water vapor. The predicted annual average contribution of 
insoluble mineral dust to CCN and CDNC in cloud-forming areas is up to 40 and 23.8%, 
respectively. Sensitivity tests suggest that uncertainties in dust size distribution and water 
adsorption parameters modulate the contribution of mineral dust to CDNC by 23 and 
56%, respectively. Coating of dust by hygroscopic salts during the atmospheric aging 
causes a twofold enhancement of the dust contribution to CCN; the aged dust, however, 
can substantially deplete in-cloud supersaturation during the initial stages of cloud 
formation and can eventually reduce CDNC. Considering the hydrophilicity from 
adsorption and hygroscopicity from solute is required to comprehensively capture the 
dust-warm cloud interactions. The framework presented here addresses this need 
and can be easily integrated in atmospheric models. 

Citation: Karydis, V. A., P. Kumar, D. Barahona, I. N. Sokolik, and A. Nenes (2011), On the effect of dust particles on global 
cloud condensation nuclei and cloud droplet number, J. Geophys. Res., 116, D23204, doi: 10.1 029/20 1 IJDO 16283. 


1. Introduction 

[ 2 ] Mineral dust is a major component of particulate matter 
in the atmosphere, accounting for more than 50% of the 
global aerosol load [Grini et al, 2005; Zender and Kwon, 
2005]. Mineral aerosol may be composed of iron oxides 
(e.g., hematite, goethite), carbonates (e.g., calcite, dolomite), 
quartz, and clays (e.g., kaolinite, illite, and montmorillonite) 
[Chou et al, 2008; Coz et al, 2009; Lafon et al, 2006; 
Twohy et al, 2009]. Differences in parent soils, and emis- 
sion and transport processes cause substantial changes 
in size distribution, composition, and morphology of dust 
particles [Jeong and Sokolik, 2007; Sokolik et al, 2001]. 
The long-range transport of dust particles can influence the 
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composition and dynamic state of the atmosphere thousands 
of kilometers downwind of their source region [e.g., Kallos 
et al, 2007; Prospero et al, 2001]. Under favorable condi- 
tions, dust particles originating from Northern and Central 
Africa may be elevated and travel across the Atlantic Ocean 
toward the Caribbean [Chiapello et al, 2005; Kallos et al, 
2006; Karyampudi and Carlson, 1988; Karyampudi et al, 
1999] or across the Mediterranean toward Europe, affect- 
ing both air quality and meteorology in Southern Europe 
[Mitsakou et al, 2008; Querol et al, 2009]. Dust from the 
Gobi and Taklamakan deserts often crosses the Pacific and 
reaches the west coast of the Americas [Meskhidze et al, 
2005a; Meskhidze et al, 2003; Sassen et al, 2002]. 

[ 3 ] Absorption of solar radiation by dust can modify the 
atmospheric thermodynamic structure, leading to suppression 
(or enhancement) of precipitation depending on cloud types 
and atmospheric conditions [Yin and Chen, 2007]. Moreover, 
dust particles are efficient ice nuclei (IN) and contribute to 
the formation of ice particles in cirrus and mixed-phase 
clouds [DeMott et al, 2003; Field et al, 2006; Teller and 
Levin, 2006; Barahona et al, 2010a]. Dust may also inter- 
act with sea salt, anthropogenic pollutants, and secondary 
organic aerosol, forming particles that consist of a “core” of 
insoluble mineral dust with coatings of soluble material 
[Gibson et al, 2007; Levin et al, 2005; Seisel et al, 2005]. 
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Dust particles with a soluble coating are typically very effi- 
cient cloud condensation nuclei (CCN), often maintaining 
their activity as IN [Levin et al., 2005]. According to Kelly 
et al. [2007], small amounts of highly soluble components 
greatly enhance the ability of fine dust particles to serve as 
CCN. Hoose et al. [2008a] have shown that, depending on 
the dust mineralogical composition, coating with soluble 
material from anthropogenic sources can lead to quasi- 
deactivation of natural dust ice nuclei. Moreover, due to their 
large size, dust particles can act as giant CCN that can form 
efficient collector drops and initiate the onset of drizzle and 
precipitation [Feingold et al., 1999; Levin and Cotton, 
2009]. The same particles can strongly compete for water 
vapor in the nucleation zone of cloudy updrafts, potentially 
reducing supersaturation and cloud droplet formation 
[Barahona et al, 2010b]. According to Lee et al. [2009], 
CCN concentration decreases in dusty regions up to 10-20% 
because dust competes for condensable H2SO4, reducing the 
condensational growth of ultrafine mode particles to CCN 
sizes. On the other hand, below 0.1% supersaturation, CCN 
concentrations increase significantly in dusty regions due to 
the presence of coarse dust particles. Manktelow et al. [2010] 
found that CCN-sized particles were depleted by an average 
of less than 5% during an extreme dust event over E. China, 
Korea, and Japan, because decreases in fine aerosols were 
mostly compensated for by increases in coarse particles. 

[4] Kohler theory (KT), which is currently used to describe 
droplet formation in atmospheric models, assumes that the 
CCN activity of aerosols depends solely on their curvature 
and the fraction of soluble material. There is emerging evi- 
dence though that insoluble species (like freshly emitted 
mineral dust) can also act as good CCN even if they lack 
appreciable amounts of soluble material. Koehler et al. 
[2009] and Herich et al. [2009] measured CCN activation 
of two types of regional dust samples (Northern Africa and 
Arizona Test Dust) and several clays (kaolinite, illite, and 
montmorillonite) at atmospherically relevant supersaturations. 
The high CCN activity was attributed to soluble ions present 
in the material. Further analysis of these data by Kumar et al. 
[2009a] showed the importance of including water adsorp- 
tion effects in describing the hygroscopic and CCN behavior 
of mineral aerosol. This was later confirmed by Kumar et al. 
[2011a, 2011b] who studied the CCN activity of dry- and 
wet-generated clays and mineral dusts representative of major 
regional dust sources (Northern Africa, East Asia/China, and 
Northern America). The observed hygroscopicity could not 
be attributed to the soluble ions present, but rather to the 
strong water vapor adsorption on the particle surface. Lathem 
et al. [20 1 1 ] also demonstrated the importance of adsorption 
for describing the water uptake properties of volcanic ash. 
Finally, Kumar et al. [201 lb] presented a “unified dust acti- 
vation framework” to treat the activation of dust internally 
mixed with soluble salts by considering the effects of 
adsorption and absorption on water activity. 

[5] The dependence of CCN activity on dry diameter dif- 
fers between particles for which soluble material dominates 
their composition and those which lack any appreciable 
amounts of solute. This difference needs to be accounted for 
in calculation of CDNC when the CCN population contains 
both particle types. Kumar et al. [2009b] developed a cloud 
droplet formation parameterization to address this need, 
which accounts for the effect of adsorption activation by 


assuming the CCN constitute an external mixture of soluble 
aerosol (that follow KT) and completely insoluble aerosol 
(that follow FHH adsorption activation theory, FHH-AT). 
Kumar et al. [2011a] also developed a framework to treat 
dust particles with substantial amounts of soluble material 
where the CCN activity is influenced by both adsorption and 
solute. 

[6] To date, global models and studies that include the 
dust impacts on cloud droplet number assume that dust acts 
as CCN only when soluble material exists in the particles 
(and controls its hygroscopicity) [e.g., Hoose et al, 2008b; 
Kelly et al, 2007; Lee et al, 2009; Manktelow et al, 2010; 
Pringle et al, 2010a, 2010b; Solomos et al, 201 1]. Here we 
assess for the first time the impact of water adsorption on the 
CCN activity of insoluble (dust) aerosol and the resulting 
contribution to global CCN and cloud droplet number. 
Simulations are carried out with the NASA Global Modeling 
Initiative (GMI) Chemical Transport Model using offline 
wind fields and an online aerosol simulation coupled with 
the Kumar et al. [2009b] activation parameterization. The 
sensitivity of the proposed modeling framework to the size 
distribution and the adsorption parameters controlling the 
dust CCN activity is investigated. Finally, the sensitivity of 
our simulations to the deposition of solute on dust particles is 
addressed by extending the Kumar et al [2009b] parame- 
terization that enables us to include the “unified” adsorption- 
absorption activation theory of Kumar et al. [2011b] for 
describing the CCN activity of aged dust particles. 

2. Model Description 

2.1. NASA Global Modeling Initiative (GMI) 

[7] The NASA GMI (https://gmi.gsfc.nasa.gov/gmi.html) 
is a modular 3-D chemistry and transport model (CTM), 
with the ability to carry out multiyear simulations for impact 
assessment studies. The detailed description of the frame- 
work is given by Rotman et al. [200 1 ] and Considine et al. 
[2005]. The GMI aerosol model used in this study was 
contributed by Liu et al [2005] and coupled to the GMI- 
CTM advection core. The aerosol module includes primary 
emissions, chemical production of sulfate in clear air and 
the in-cloud aqueous phase, gravitational sedimentation, dry 
deposition, wet scavenging in and below clouds, and hygro- 
scopic growth. The model time step for chemistry is one 
hour. Model inputs include emissions of SO2 (fossil fuel and 
natural), DMS, H2O2, black carbon (from biomass burning 
and fossil fuel), organic carbon (fossil fuel, biomass burning, 
and natural), mineral dust (four size bins), and sea salt (four 
size bins), which are provided by Liu et al. [2005]. The 2% 
(by mole) of fossil-fuel sulfur emission is assumed to occur 
as primary sulfate aerosol to account for the rapid conversion 
of SO2 to S04^ in combustion plumes. SO2 emissions also 
include sources from sporadic and continuously emitting 
volcanoes (averaged over a 25-year time period). Biofuel 
emissions are categorized as fossil fuels. 

[s] Dust aerosol fluxes at every 6 h were generated with 
the approach of Ginoux et al. [2001], using the NASA 
Goddard Data Assimilation Office (DAO) meteorological 
fields for year 1997. To account for a systematic overesti- 
mation in calculated AOD over deserts, Liu et al. [2005] 
reduced the submicron dust (0.1-1.25 ^m bin) emission 
flux by roughly twofold. With this modified emissions 
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Table 1. Size Distribution Parameters Applied to Fossil Fuel, Biomass Burning, Marine, and Mineral Dust Aerosols 


Aerosol Type 

Aerosol 

Components 

Density 
(g cm^^) 

Median 

Diameter 

(pm) 

Geometric 

Standard 

Deviation 

Number 

Fraction 

Reference 

Fossil fuel 

sulfate 

1.77 

0.1 

1.9 

1 

Chuang et al. [1997]; Radke et al. [1988] 


organic carbon 

1.2 






black carbon 

1.5 





Biomass burning 

organic carbon 

1.2 

0.16 

1.65 

1 

Anderson et al. [1996] 


black carbon 

1.5 





Marine 

natural sulfate 

1.77 

0.018 

1.4 

0.81 

Lance et al. [2004] 


sea salt 

2.2 

0.075 

1.6 

0.18 





0.62 

2.7 

0.01 


Mineral dust (base case) 

mineral dust 

2.6 

0.16 

2.1 

0.93 

d’Almeida et al. [1987] 




1.4 

1.9 

0.07 





10 

1.6 

3 X 10“*^ 


Mineral dust (sensitivity case) 

mineral dust 

2.6 

0.18 

1.8 

0.76 

Chou et al. [2008] 




0.7 

1.4 

0.13 





1.5 

2.0 

0.11 



scheme, Liu et al [2005] predicted a total dust flux (0. 1-20 ^m) 
of 1684 Tg yr^^, which is on the lower end of predietions 
from ten different models (820-5102 Tg yr^^). Liu et al. 
[2005] carried out an evaluation of the aerosol simulation 
against in situ observations from 15 oceanic sites. In regions 
of the southern hemisphere, which are influenced by dust 
emissions from Patagonian (Palmer Antarctic) and Australian 
(Cape Grim, Norfolk, New Caledonia) deserts, Liu et al. 
[2005] found that the model reproduces the seasonal cycle 
of dust concentration, although the dust concentration was 
underpredicted during the austral summer. In the Northern 
Hemisphere, the model reproduced the concentrations and 
seasonal cycles of the African dust sampled at Izana Tenerife, 
Barbados, and Bermuda; the model overpredicted, however, 
the concentration in August at Izana and underpredicted the 
concentration in July and August at Bermuda by a factor of 
two. Further details on the GMI framework can be found in 
Barahona et al. [2010a, 2011]. 

2.2. Aerosol Distribution 

[ 9 ] The concentration of particles that can experience 
hygroscopic growth is given as an input from the aerosol 
module to the cloud droplet formation parameterization and 
is distributed in four aerosol types: fossil fuel (sulfate, 
organic mass, and black carbon), biomass burning (organic 
mass and black carbon), marine (natural sulfate and sea salt), 
and mineral dust. Fossil fuel, biomass burning, and marine 
aerosols are assumed to follow KT for CCN activation, 
whereas mineral dust is assumed to be insoluble and follow 
FHH-AT. Particles within each aerosol type are internally 
mixed and assumed to follow a prescribed number size dis- 
tribution function (Table 1). Fossil fuel aerosols are repre- 
sented by a size distribution given by Chuang et al. [1997] 
and Radke et al. [1988]. The particle size distributions for 
biomass burning aerosols are based on measurements by 
Anderson et al. [1996]. Marine aerosols are represented by a 
size distribution proposed by Lance et al. [2004] and mineral 
dust particles are assumed to follow the size distribution by 
d’Almeida et al. [1987]. The number concentration of each 
aerosol type (e.g., fossil fuel) is given by A = mtJ{Vptot), 
where nitot is the total mass concentration (across size bins) 
of the aerosol type, V is the volume of a particle with dry 
diameter equal to dg, and p,ot is the density of the aerosol 


type given by ptot = m,ot Y. /ft)- Here m,- is the mass 

/ i=i 

concentration of each component within the aerosol type 
(e.g., sulfate, OC, BC), p,- is its density, and n is the number 
of aerosol types. 

2.3. Cloud Droplet Formation Parameterization 

[ 10 ] Calculation of CDNC for particles following FHH-AT 
is performed in two conceptual steps, one involving the 
determination of the “CCN spectrum” (i.e., the number of 
CCN that can activate to form droplets at a certain level of 
supersaturation), and another one determining the maximum 
supersaturation, that develops in the ascending cloudy 
parcels used to represent droplet formation in the general 
circulation model (GCM). The CDNC is then just the value 
of the CCN spectrum at 5max- 

[ 11 ] The “CCN spectrum,” F‘\s), is computed following 
Kumar et al. [2009b] and assumes that particles can be 
described either by KT or FHH-AT. F‘\s) for an external 
mixture of lognormal particle size distributions is given by: 


P(s) 




X'/l ln((T;) 


( 1 ) 


where s is the level of water vapor supersaturation, Sg i is the 
critical supersaturation of the particle with a diameter equal 
to the geometric mean diameter of the mode i, cr, is the 
geometric standard deviation for mode i, and x is an expo- 
nent that depends on the activation theory used. For modes 
following Kohler theory, x = -3/2 [Fountoukis and Nenes, 
2005], while for particles following FHH theory, x depends 
on ApHH and Bphh [Kumar et al, 2009b]. The adsorption 
parameter Bp^^ expresses the long range interactions of 
adsorbed water layers with the dust surface. The smaller the 
value of BpHH, the greater the distance from the dust surface 
for which attractive forces are present. Ap^^ primarily 
represents the interaction between the first water monolayer 
and the dust surface. Appipi and Bp^n are compound-specific 
and determined experimentally, with Ap^H typically ranging 
from 0. 1-3.0, and Bp^^ ranging from 0. 5-3.0 [Sorjamaa 
and Laaksonen, 2007]. Kumar et al. [2011b] tested a wide 
range of fresh unprocessed regional dust samples and 
minerals and found that one set of the FHH parameters 
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i^FHH = 2.25 ± 0.75, BpHH ^ 1.20 ± 0.10) can adequately 
reproduce the measured CCN activity for all dust types 
considered. 

[ 12 ] The maximum supersaturation, 5max, in the ascending 
parcel is ealculated from an equation that expresses the 
supersaturation tendency in cloudy air parcels, which at the 
point of maximum supersaturation becomes \Barahona and 
Nenes, 2007; Nenes and Seinfeld, 2003], 

2aV 

- Gjmax7(0,5max) = 0 (2) 

where V is the updraft velocity, is the density of water, 
and a, 7 , G are parameters defined in Barahona and Nenes 
[2007]. Once 5max is determined by numerically solving 
equation (2), from equation (1). 

[ 13 ] In equation (2), /(0,5max) is known as the “eonden- 
sation integral” [Barahona and Nenes, 2007; Barahona 
et al, 2010b; Kumar et al., 2009b] and it expresses the 
condensational depletion of supersaturation upon the grow- 
ing droplets at the point of 5max in the cloud updraft. It is 
expressed as the sum of three terms [Kumar et al., 2009b], 

7(0Amax) — S max ) + Ifhh[0,s max ) + Igccn (0, ^ max ) (3) 


The first term in the right hand side of equation (3) gives 
the eontribution from particles that follow Kohler theory, 
Ihe second one from particles that follow the 
FHH theory, and the third one accounts for 

the effect of the very large and giant CCN (neglected by the 
other two terms). Using the population splitting approach 
of Nenes and Seinfeld [2003], 7^(0, 5max) is calculated as: 

*Smax) ^ ^K.li^D^part) 3“ dp2i.^parTf^maf)'> where dp i{Q,Spart^ 

and iK.iiSpart, Smax) are presented by Fountoukis and Nenes 
[2005] or Kumar et al. [2009b]. The partition supersatura- 
tion, Spart, separates two CCN populations, one (expressed by 
kK, 2 {Spart,Smaf)) for which droplets experience negligible 
growth beyond the critical diameter ( 5 ^ = ^max), and another 
one (expressed by lK.\{^,Spart)) for which droplet growth is 
much larger than the critical diameter {Sc ^ •Smax)- In 
equation (3), lFHH{^,Smaf) represents the contribution of 
FHH particles to the condensation integral [Kumar et al., 
2009b]. 

[ 14 ] If dust particles contain a substantial fraction of soluble 
material, the critical supersaturation is determined using the 
“unified dust activation framework” of Kumar et al. [201 lb], 
which involves determining the maximum of the relevant 
equilibrium curve describing an aerosol particle (consisting 
of insoluble core with a soluble coating) in equilibrium with 
the surrounding water vapor: 


4(tM„ 
RT p„.Dp 


EsD. 


Dj, 


dry 

F 

^i^dry 


— Ap 


Dp - eJ^D, 

2D,, 


dry 


~Bfhh 


(4) 


where Dp is the wet equilibrium diameter, D^i^y is the particle 
dry diameter, is the insoluble volume fraction, is the 
soluble volume fraction, and k is the hygroscopicity of the 
soluble fraction. The exponent x in the case of aged dust 
particles lies somewhere between the KT and FHH-AT 
limits, and is determined (as described by Kumar et al. 
[2011b]) by performing a power law fit between Sc and 


D,iry determined from equation (4). Given that dust is 
considered to be a separate population from other aerosol 
types, application of the unified dust activation framework 
in the Kumar et al. [2009b] parameterization is carried out 
by replacing the x, s^ normally computed by FHH-AT 
with those of the unified framework. 

[ 15 ] Finally, /Gccv(0>5max) accounts for the effect of the 
very large and giant CCN on the condensation rate and is 
presented by Barahona et al. [2010b] as: 


foCCNi^T ^ ^ I 2 Dg,i 

■ erfc 


/ V \ 2 


1 , X 

In -A V—f=\nGi 


x\[2 Ino-,' Dg i x[2 


(5) 


where Dp^^^ is the wet equilibrium diameter of the smallest 
particle that activates (i.e., for which Sc = ^max) and Dg i is the 
geometrie mean diameter of mode i. For giant nuclei following 
Kohler theory, Dp^^ = 2A/(3^/3spart), where A = 4aMJRTp„, 
and represents the boundary between partieles that experi- 
ence significant growth after activation and those that are 
strongly kinetically limited [Barahona et al., 2010b; Nenes 
and Seinfeld, 2003]. Giant nuclei following FHH theory do 
not grow significantly after activation and therefore are 
kinetically limited by the inertial mechanism [Nenes et al, 
2001]; indeed, DJD^iry < 2 for most atmospherically rele- 
vant combinations of Afhh and Bfhh [Kumar et al., 2009b], 
where is the critical wet droplet diameter. Dp^^ is calcu- 
lated at Spart, given by the solution of equation (4). 

2.4. Model Application 

[ 16 ] The model results presented here are based on 1 yr 
simulations and the first month of each simulation has been 
excluded in order to limit the effeet the initial conditions 
have on the results (simulations with longer startup periods 
did not substantially affect the simulation results). The 
meteorological fields used in the simulations were taken 
from the Goddard Institute for Space Studies version II 
(GISS II) GCM [Koch and Rind, 1998; Rind and Lerner, 
1996], which includes a slab (Q-flux) ocean model to repre- 
sent the ocean-atmospheric coupling. The data set spans 
over 1 yr and represents the period from March 1997 to 
February 1998. The horizontal resolution is 4° latitude by 
5° longitude. The vertical resolution is 23 vertical layers 
(from surface to 0.017 hPa). Meteorological information 
was updated at 3 h intervals and archived for the year of 
simulation. 

[ 17 ] Parameters used by the cloud droplet formation 
parameterization in this study include the aerosol size dis- 
tributions (section 2.1), an effective water vapor uptake 
coefficient of 0.06 [Fountoukis et al, 2007], the updraft 
velocity representative of typical stratocumulus clouds and 
constrained using observations, F= 0.3 m s ' over land, and 
F = 0.15 m s^' over ocean [Chuang et al, 2000; Guibert 
et al, 2003; Meskhidze et al, 2005b], and “base case” 
FHH adsorption parameters, Afhh ^ 2.25 and Bfhh =1.2 
[Kumar et al, 2011b]. Updraft velocity varies significantly 
within a GCM grid cell and between clouds; however, 
grid-average CDNC for stratiform clouds can be calculated 
by computing CDNC at a “characteristic” vertical velocity 
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Figure 1. Predicted annual (top) mean cloud droplet and (bottom) cloud condensation nuclei (CCN) at 
0.2% supersaturation concentrations (cm ^) for the lowest cloud-forming level (960 mbar). White repre- 
sents cloud-free areas. 


[Morales and Nenes, 2010], which is very close to the 
average updraft velocity. Meskhidze et al. [2005b] and 
Fountoukis et al. [2007] showed that using the average 
updraft gives optimal CDNC closure in cumulus and strato- 
cumulus clouds sampled in a variety of field campaigns. 

3. Results 

3.1. Overview of Model Predictions 

[is] The predicted annual mean CCN at 0.2% supersatu- 
ration and cloud droplet number concentrations for the 
lowest cloud-forming level (at 960 mbar) are shown in 
Figure 1. CDNC is referred to the number concentration of 
droplets nucleated in-cloud and represents an upper limit 
since the model does not account for droplet depletion by 
collision, coalescence and collection. As expected, higher 
CDNC concentrations (up to 1000 cm^^) are predieted over 
the midlatitudes of the Northern Hemisphere (i.e., over 
China, Europe, and eastern U.S.; Figures la and 2), consis- 
tent with the high concentration levels of CCN associated 


with industrialized regions (Figure lb). Over the eon- 
tinents of the Southern Hemisphere, large CDNC values 
(up to 600 cm^^) occur over regions affected by biomass 
burning in South America and Africa. The lower continental 
CDNC values are predicted over the arid areas of the Sahara 
and Gobi deserts (up to 300 cm~^). Over oceans, CDNC is 
increased up to 300 cm~^ by continental aerosol transported 
by the trade winds off the subtropical west eoasts of Africa 
and America, and westerlies in midlatitude east eoasts of 
North America and Asia (Figure 1) [Minnis et al, 1992; 
Prospero et al., 1983]. In contrast to the Northern Hemi- 
sphere, lower cloud droplet number concentrations are found 
over the cleaner remote oceans of the Southern Hemisphere 
and in the polar regions (Figure la, lb). In these regions, 
CDNC is predicted to be up to 40 em Over the Southern 
Ocean, larger CDNC values are predicted (up to 80 ciW^) 
due to the enhaneed production of sulfate from dimethyl 
sulfide (DMS) oxidation, and the high concentrations of sea 
salt, which is lofted by winds associated with the Antarctic 
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Figure 2. Predicted zonal annual mean cloud droplet num- 
ber concentration (cm^^). White areas correspond to non- 
cloud-forming regions. 



circumpolar vortex that flows eastward around Antarctica. 
Predicted CDNC decreases with altitude over the polluted 
areas of the Northern Hemisphere caused by a decrease in 
aerosol concentration (Figure 2). Over the tropics, however, 
CDNC increases until 865 mbar, influenced by long-range 


transport of polluted air masses, and then decreases. The 
global mean CDNC is predicted to be 169 cm^^, 177 cm^^, 
188 cm”^, 181 cm”^, and 128 cm”^ at 960, 920, 865, 780, and 
650 mbar, respectively. These values are within the limits of 
simulated CDNC published in the literature. Merikanto et al. 
[2010] obtained a yearly global average CDNC of 211- 
240 cm~^ from model rans at 300-1000 m above ground 
level. Chen etal. [2010] reported a CDNC value of 189 cm^^ 
at 972 mbar and 122 ciW^ at 850 mbar. Penner et al. [2006] 
found the CDNC ranging from 119 cm~^ to 159 cm^^ near 
850 mbar. Barahona et al. [2011] found the global mean 
CDNC ranging from 96 cm^^ to 103 cm^^ at 850 mbar 
depending on the climatological meteorological fields used. 
Finally, Leibensperger et al. [2011] report a global mean 
CDNC of 83 cm~^ at 850 mbar. 

3.2. Model Evaluation 

[ 19 ] The predicted cloud droplet number concentration is 
compared against observational data from continental, pol- 
luted marine, and clean marine regions around the world 
(Figure 3). The predictions reflect aerosol distributions in an 
average climate state; given that the observations span over a 
decade, we expect a spatial and seasonal correspondence with 
the predictions. The summary of this comparison is depicted 
in Table 2 and Figure 4. Over clean marine regions. 



Figure 3. Location of observational data used for the evaluation of model predictions. Continental, pol- 
luted marine, and clean marine regions are represented by red, green, and blue symbols, respectively. 
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Figure 4. Comparison of simulated and observed global cloud droplet number concentrations (cm^^). 
Open symbols correspond to dust affected regions. Also shown are the 1:1, ±50% lines, and the probabil- 
ity distribution of the ratio of the simulated CDNC to the observed CDNC (inset plot). 


concentrations of 27 cm”^ (South Pacific) to 108 cm~^ 
(Eastern Pacific) are predicted; this agrees with retrieved 
CDNC values mostly to within 50%, ranging from 40 cm ^ 
(South Pacific [Bennartz, 2007]) to 279 cm^^ (Eastern Pacific 
[Straub et al, 2007; vanZanten et al, 2005]). Larger under- 
estimation is found over the remote area west of Australia 
(41 cm'^ versus 107 cm~^ reported by Rausch et al. [2010]). 
Moreover, over the Western Arctic Ocean, the predicted 
CDNC range from 158 cm~^ to 200 cm ^ . This agrees well 
with available observations [Herman and Curry, 1984; 
Hobbs and Rangno, 1998], but is considerably overpredicted 
over the North Pacific and Arctic Ocean [Gultepe and Isaac, 
2002; Gultepe et al., 2004; Luo et al., 2008]. 

[ 20 ] The predicted CDNC increases to 88-580 cm'^ over 
polluted marine regions close to coasts. The highest observed 
CDNC value is over the Yellow Sea (eastern coast of China). 
In particular, Zhao et al. [2006] reported measured values 
ranging Irom 30 cm”^ to 1000 cm~^. The predicted CNDC is 
226-263 cm"^ which is relatively high compared to satellite 
retrievals (100-250 cm~^ reported by Bennartz [2007] and 
Rausch et al. [2010]). The model predictions are in better 
agreement with these studies over North Ameriea coast 
(124 cm^^ versus 96-300 cm^^). South America coast 
(210 cm~^ versus 77-300 cm^^), and Africa (1 13-161 cm”^ 
versus 95-300 cm^^). Over Atlantic Ocean, the CDNC pre- 
dictions range between 88 cm~^ (eastern north Atlantic) and 
240 cm~^ (Santa Maria, Azores). This is well within the 
observed values over that region (65-300 cm^^ [Brenguier 
et al, 2003; Chuang et al, 2000; Dong et al, 1997; 
Harshvardhan et al, 2002; Schuller et al, 2005]). 

[ 21 ] The predicted CDNC over continental regions ranges 
between 201 cm”^ (Barrow, AK) and 956 cm^^ (Koblenz, 
Germany) compared to observed values of up to 2000 cm^^ 
(Cumbria, North England [Bower et al, 1999]). In Asia, 
predicted values range Irom 577 cm^^ over Western China 
(observed 30-700 cm”^ [Zhao et al, 2006]) to 772 cm~^ 


over Beijing (observed 30-1000 cm~^ [Zhao et al, 2006]). 
In Europe, the model predicts a lower concentration 
(396 cm ^) over North England, where there are the highest 
observations (100-2000 cm^^ [Bower et al, 1999]). The 
highest CDNC are predicted over Koblenz, Geraiany (804- 
956 cm ^) and are in good agreement with observations 
(675-965 cm~^ [Lehmann et al, 2009]). Over Finland, the 
model tends to overpredict CDNC, with values ranging from 
550 cm^^ to 899 cm^^ versus 30 cm ^ to 610 cm ^ [Komppida 
et al, 2005; Lihavainen et al, 2008; Portin et al, 2009]. 
Finally, over North America, the predieted CDNC range from 
201 cm~^ over Alaska (observed of 54—222 cm~^ [Dong and 
Mace, 2003]) to 896 em“^ over Oklahoma (observed of 
100-650 cm~^ [Dong et al, 2005, 1997, 2000; lacobellis 
and Somerville, 2006]). The observations range from 
54 em^^ (Barrow, AK [Dong and Mace, 2003]) to 1300 cm^^ 
(Cleveland, OH, and Detroit, MI [Fountoukis et al, 2007]). 

[ 22 ] Overall, the model is in reasonably good agreement 
with observations over both clean and polluted marine 
regions as well as over regions that are relatively close to 
dust sources (e.g., Canary Islands). More precisely, 67% of 
the CDNC predictions over those regions diverge less than 
50% from available measurements (Figure 4). The model 
tends to overprediet the CDNC compared to observed data 
from continental regions (Figure 4). Considering the influ- 
ence of droplet collision and coalescence processes may, in 
part, reduce prediction biases. Given however that dust is not 
expected to be a significant contributor to CCN number in 
these regions, we do not anticipate the overprediction bias to 
affect conclusions. 

3.3. Effect of Dust on CCN and CDNC 

[ 23 ] The contribution of insoluble particles to cloud 
droplet formation is potentially most important in areas close 
to mineral dust sources. These areas may not always be 
associated with extensive cloud cover; to account for this. 
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Mean 1 % Max 40% 



0 4 8 15 30 50 


2 6 10 20 40 

Figure 5. Predicted annual mean contribution (%) of min- 
eral dust to particle number concentrations for the lowest 
cloud-forming level (960 mbar). White represents areas that 
are free of wann clouds throughout the year. 


[ 25 ] The contribution of dust to CNDC is computed 
directly from the parameterization output, as the activation 
fraction for each aerosol mode is known. The predicted Srrmx 
in clouds forming in these dusty areas remains low (~0.1% 
on average), hence mineral dust contributes up to 23.8% of 
the CDNC (Figure 7). Moreover, allowing dust to act as 
CCN resulted in a relatively significant decrease of 5„ax over 
Western Saharan (30%) and Patagonian (10%) deserts. Over 
the remaining arid areas (e.g., Gobi and Australian deserts) 
the decrease in s^ax is less important (less than 5%). Mineral 
dust also has a small contribution to the annual average 
predicted CDNC across the Atlantic Ocean, as far as the 
Caribbean Sea (up to 10% near the North African coast). 
This contribution can be significantly larger during dust 
episodes. The peak contribution of mineral dust to CDNC 
across the Atlantic Ocean is predicted during the month of 
Febmary (15%). The greatest seasonal variability of the con- 
tribution of mineral dust to CDNC is predicted downstream 
of the Patagonia Desert, where it varies from 2% during May 
to 40% during February. 


CDNC are calculated only in cells for which warm clouds 
are present (i.e., cloud water is present, and cloud top tem- 
perature exceeds 263K). Mineral dust effects on CDNC 
depend on the predicted contribution of mineral dust to total 
aerosol number concentration, as the available aerosol bud- 
get is the main factor driving the CCN population. Close to 
deserts, the dominant aerosol constituent (80-99% by mass) 
is mineral dust with the rest of the aerosols consisting of 
soluble materials, such as sulfate, organic carbon, and sea salt 
from long-range transport. However, these soluble aerosols 
are small in size compared to mineral dust and therefore can 
contribute substantially to aerosol number and CCN con- 
centrations. The predicted contribution of mineral dust to 
particle number concentrations is presented in Figure 5. The 
contribution of dust to the total aerosol number concentration 
is small on a global scale (~1%), but regionally can be large, 
especially during dust episodes and their long-range transport 
(Figure 5). The largest contribution (up to 40%) is seen over 
North Africa and Asia (Arabian Peninsula and Gobi Desert) 
while mineral dust from South America (Patagonian Desert) 
and Australian deserts also contributes to the total aerosol 
number concentration (10-15%). These contributions repre- 
sent an average climatology but will vary significantly, both 
seasonally and interannually. 

[ 24 ] To investigate the potential importance of dust on 
cloud droplet formation, the predicted contribution of dust to 
total CCN is shown in Figure 6 for moderate (0.2%) and 
high (0.4%) supersaturations. The contribution of mineral 
dust to the CCN concentration decreases with decreasing 
supersaturation as insoluble aerosols are less hygroscopic 
and they need higher supersaturations in order to activate 
compared to soluble particles. Nevertheless, dust particles 
are coarser than soluble aerosols and therefore, if coated 
with hygroscopic material, they can potentially activate at 
lower supersaturations than smaller anthropogenic particles. 
At 0.2% supersaturation, the contribution of mineral dust to 
the total CCN concentration is up to 30% while at 0.4% 
supersaturation the contribution increases to 35%. In both 
cases, the maximum predicted contribution occurred on the 
eastern part of the Arabian Peninsula. 


4. Sensitivity Analyses 

[ 26 ] Dust-relevant inputs in the CDNC parameterization, 
other than dust concentration, include Apfm, Bp^H and the 


Mean 0.4% 


Max 30.1% 



-180 -120 -60 0 60 120 180 
0 4 8 15 30 50 


Mean 0.6% Max 35.1% 



Figure 6. Predicted aimual mean contribution (%) of mineral 
dust to total CCN for (top) 0.2% and (bottom) 0.4% super- 
saturation for the lowest cloud-forming level (960 mbar). 
White represents areas that are free of wann clouds through- 
out the year. 
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Mean 0.3% Max 23.8% 



-180 -120 -60 0 60 120 180 


0 2 4 6 8 10 20 

1 3 5 7 9 15 


Mean 3.7% Max 88.8% 



I 


5 15 25 40 60 100 

Figure 7. Predicted annual mean mineral dust contribution 
(%) to cloud droplet number concentration for the lowest 
cloud-forming level (960 mbar) by using (top) present-day 
and (bottom) preindustrial anthropogenic emissions. White 
represents areas that are free of warm clouds throughout 
the year. 


size distribution of dust particles. To investigate the impact 
of uncertainty in these inputs on the predicted effect of 
mineral dust on CDNC, five sensitivity simulations were 
conducted. First, we carry out simulations using an alternate 
size distribution (Table 1), measured over Niger during 
January 2006 [Chou et al, 2008]. Second, we vary the 
sensitivity to Ap^H and Bp^H parameters using limits 
reported by Kumar et al. [2011b]. In particular, the set of 
parameters used in these simulations are: Apupi = 2.25, 
BpHH =1-1; Apiiii - 2.25, Bpi^fj = 1.3; Apfjf/ ~ 1.5, Bpi^p/ = 
1.2; and Ap^n = 3.0, Bp^n = T2. Finally, two additional 
sensitivity simulations were conducted to test the effect of 
anthropogenic emissions and atmospheric aging on mineral 
dust contribution to CDNC. 

4.1. Sensitivity of CDNC and Total Aerosol Number 
to Size Distribution of Dust 

[ 27 ] The size distribution used for the sensitivity test tends 
to distribute more mineral dust mass onto larger size parti- 
cles than in the base case scenario. As a result, the number 
concentration of mineral dust is lower; compared to the base 
case simulation, the contribution of insoluble aerosols to 
annual mean total aerosol number concentration is found to 


be 45% lower (not shown). As expected, this significant 
change in the number of available insoluble particles reduces 
the contribution of dust particles to CDNC. The fraction of 
the droplets coming from the activation of insoluble particles 
is on average 23% lower than the base case simulation in 
areas with high concentration of dust (Table 3). Over remote 
areas of the Atlantic Ocean, the contribution of mineral dust 
to CDNC is 40% lower than the base case simulation (not 
shown). The base case simulation predicted that mineral dust 
can contribute to cloud droplet formation over a wide area 
that extended thousands of miles from its sources (e.g., the 
tropical Atlantic Ocean), in contrast to the sensitivity case 
simulation where the mineral dust effect on CDNC is pre- 
dicted in a relatively limited area close to its sources (i.e., 
deserts). 

4.2. Sensitivity of CDNC to and tS/r//// Parameters 

[ 28 ] Lower values of the Bp^^ parameter corresponds to a 
more hydrophilic particle and largely determines the exis- 
tence and values of the particle critical supersaturation, 
[Kumar et al., 2009b]. Appipi also affects the particle CCN 
activity, but to a lesser extent than Bp^n- Kumar et al. [2011b] 
showed that one set of the FHH parameters (Apun ^ 2.25 ± 
0.75, Bpipj = 1.20 ± 0.10) encompasses the range of hygro- 
scopicity observed for a wide range of dry-generated dust. 
The upper and lower bounds of this set were used as sensi- 
tivity tests in order to determine the effect of the level of 
hygroscopicity of dust on CDNC. When Bp^m = 1.3, the 
maximum contribution of insoluble particles to the predicted 
CDNC decreased to 15.1% as the mineral dust become less 
hygroscopic compared to the base case simulation (Table 3). 
Using the lower bound Bp^^ = 1.1, almost doubled the 
maximum contribution of dust on CDNC to 37.2% (Table 3). 
The set of the FHH parameters measured by Kumar et al. 
[2011a] represents the hygroscopicity of freshly emitted dust. 
Nevertheless, during the aging of mineral dust, its hygro- 
scopicity increases as it is coated with soluble materials (e.g., 
sea salt or sulfates), and therefore the value of 1 . 1 should not 
considered as an upper limit of hygroscopicity, as a freshly 
emitted dust particle with diameter 100 nm and Bppm = 
1 . 1 exhibits comparable hygroscopicity to a KT CCN with a 
(NH 4 ) 2 S 04 volume fraction of 10% [Kumar et al., 2011a]. 
Finally, the results are less sensitive to changes in Ap^n, 
as the contribution of insoluble particles on the predicted 
CDNC varies from 17.3% (when Appifj = 1.5) to 27.4% 
(when Apfjfj =3.0) with the base case value being equal to 
23.8% (when Ap^n = 2.25; Table 3). 


Table 3. Annual Mean Contribution of Mineral Dust to Cloud 
Droplet Number Concentrations Predicted by the Base Case and 
the Sensitivity Case Simulations 


Simulation Scenario 

Contribution 

Base case 

23.8% 

Sensitivity to size distribution 

18.0% 

Sensitivity to Bphh (Bphh= 1-3) 

15.1% 

Sensitivity to Bphh 1-1) 

37.2% 

Sensitivity to Afhh (^fhh = 3.0) 

27.4% 

Sensitivity to Afhh (Afhh= 1.5) 

17.3% 

Sensitivity to anthropogenic 

88.9% 

emissions 
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Figure 8. Predicted annual mean cloud droplet concentra- 
tions (cm activated from (top) completely insoluble dust 
particles and (bottom) aged dust particles for the lowest 
cloud-forming level (960 mbar). 

4.3. Sensitivity of Mineral Dust Contribution to CDNC 
to Anthropogenic Emissions 

[ 29 ] Current emission control policies are designed to 
lower the anthropogenic fine particulate matter concentration 
in order to meet the air quality standards. This will enhance 
the importance of mineral dust to total aerosol budget and 
cloud droplet formation over areas affected by anthropogenic 
emissions. To estimate an upper bound of mineral dust con- 
tribution to cloud droplet formation, a simulation excluding 
the anthropogenic emissions (emissions from fossil fuels) is 
conducted (pre-industrial conditions). In that case, the model 
predicted that close to deserts, up to 88.8% of cloud droplets 
formed from the activation of insoluble particles (Figure 8). 
On global average, the predicted contribution of mineral dust 
to CDNC is 12 times higher than the base case simulation 
(3.7% instead of 0.3% in the base case). 

4.4. Sensitivity of Mineral Dust Contribution to CDNC 
to Atmospheric Aging 

[ 30 ] During atmospheric transport, fresh dust undergoes 
aging, which results in a coating of soluble material on its 
surface that augments its CCN activity. In order to address 
the first order effect of aging to dust CCN activity and 
CDNC, a sensitivity study is conducted where dust is 
assumed to be coated with 10% (in volume) ammonium 


sulfate (called hereafter “aged dusf ’) and it follows the same 
size distribution as in the base case. The “unified dust acti- 
vation framework” of Kumar et al. [2011a] (section 2.3) is 
used to calculate the CCN activity of aged dust; all other 
parameters are kept as in the “base case” simulation. Com- 
parison between CDNC coming from the activation of 
completely insoluble dust (Figure 8a) and those from aged 
dust particles (Figure 8b) suggests that in the latter case the 
contribution of dust to CDNC increased by 2.5 times on 
average. The maximum CDNC (predicted over the Gobi 
Desert) inereased from 38 cm~^ to 1 15 cm"^. In contrast with 
the completely insoluble dust particles (see section 3.3), the 
contribution of aged dust to CCN concentration increase with 
decreasing supersaturation. The increased hygroseopicity of 
aged dust together with its large dry diameter results in 
its activation at lower supersaturations than the smaller 
anthropogenie particles. At 0.4% supersaturation, the global 
annual mean contribution of mineral dust to total CCN con- 
centration is 1% (with a maxima of 63%) while at 0.2% 
supersaturation the global annual mean eontribution increa- 
ses to 1.3% (with a maxima of 75%). 

[ 31 ] Compared to freshly emitted dust, aged dust aetivates 
earlier in a cloudy updraft and grows to a larger wet size for 
a given supersaturation [Kumar et al., 2011b]; this can 
promote water uptake during the initial stages of cloud for- 
mation, eventually leading to supersaturation (and CDNC) 
depression. This is illustrated in Figure 9, whieh depiets the 
fractional change of CDNC between a simulation assuming 
no aetivation of dust and a simulation assuming activation of 
aged dust. The total CDNC is decreased by as much as 70% 
(Gobi Desert), 60% over the Arabian Peninsula, 40% over 
the Sahara Desert, 20% over Australian desert, 5% over the 
Southwestern United States, and increased by up to 20% 
over the Patagonian Desert. The CDNC also decreases over 
a wide range of areas sueh as Southern Europe (up to 5%), 
Mediterranean Sea (up to 20%), tropical Atlantic Ocean 
(up to 20%), etc. As expected, j^ax is also reduced signifi- 
cantly and the maximum decrease is observed over Arabian 
Peninsula and the Gobi Desert (up to 80%). The ^max 
decreases by up to 40% over the Sahara Desert, 30% over 


Mean -1.5% Min -67.6% Max 18.0% 
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Figure 9. Predicted annual mean fraetional change (%) of 
cloud droplet number concentration for the lowest cloud- 
forming level (960 mbar) between a simulation assuming 
no activation of mineral dust and a simulation assuming acti- 
vation of aged dust. 
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the Australian and Patagonian deserts, 5% over the South- 
western United States and Southern Europe, 20% over the 
Mediterranean Sea, and 30% over the tropical Atlantic 
Ocean. 

5. Conclusions 

[ 32 ] This study is a first attempt to assess the contribution 
of freshly emitted insoluble dust particles to global CCN and 
cloud droplet number concentrations. Simulations are car- 
ried out with the NASA Global Modeling Initiative (GMI) 
Chemical Transport Model using wind fields computed with 
the Goddard Institute for Space Studies (GISS) general cir- 
culation model. GMI simulates global distributions of fossil 
fuel, biomass burning, marine, and dust aerosols; particles 
within each aerosol type are internally mixed and assumed 
to follow a prescribe size distribution. Particles with a sig- 
nificant amount of soluble material (the CCN activity of 
which is given by Kohler theory) are assumed to be sulfate, 
organic carbon, and sea salt. Dust is assumed insoluble, 
with the CCN activity given by adsorption activation theory; 
parameters required to constrain the theory are obtained from 
activation experiments of resuspended desert soil samples 
[Kumar et al, 2011b]. Calculation of droplet number from 
the aerosol simulation is done online in GMI, using the 
Kumar et al. [2009b] parameterization, which considers 
cloud droplet formation within an ascending air parcel con- 
taining an external mixture of soluble particles (that activate 
according to Kohler theory), and insoluble, wettable particles 
(that form droplets through adsorption activation). This new 
framework is used to assess the impact of dust and adsorption 
activation on the predicted droplet number concentration. 

[ 33 ] The contribution of mineral dust to total particle and 
cloud droplet number concentration is relatively important 
in areas with high dust concentrations (e.g., deserts) with 
impacts that can extend across the tropical Atlantic Ocean. 
The predicted contribution of mineral dust to particle number 
concentrations is up to 40% over North African and Asian 
deserts and up to 15% over South American and Australian 
deserts. The effect of mineral dust on the annual average 
CDNC is up to 23.8% while the impact on the predicted 
CDNC across the tropical Atlantic Ocean is up to 10% close 
to the shore of the Western Sahara, up to 5% to the ocean, 
and less that 2% to the Caribbean Sea. Given that all the 
results in this study are expressed as annual averages, this 
contribution can be even more important during specific dust 
storm episodes. It is worth mentioning though that the model 
does not account for aerosol nucleation in the upper tropo- 
sphere and thus underestimates the aerosol number concen- 
tration. This underestimation also affects the influence of 
dust on total aerosol number concentration and consequently 
on CDNC. 

[ 34 ] Sensitivity tests show that the results are sensitive to 
the level of hygroscopicity of dust, and a 10% change in the 
BpHH parameter can lead to up to 40% change in the pre- 
dicted number of activated insoluble particles. Moreover, the 
results are also sensitive to the size distribution of dust. 
Using a size distribution of dust, measured over Niger by 
Chou et al. [2008], which distributes dust in larger particles 
than the base case distribution (by d’ Almeida [1987]), the 
predicted mineral dust particles decreased by 45% and the 


predicted cloud droplets activated from these insoluble par- 
ticles by 20%. 

[ 35 ] As current emission control policies aim to reduce the 
concentration of anthropogenic aerosol, the impact of dust 
on global aerosol budget and cloud droplet formation will 
become larger in the future. The theoretical upper bound of 
the potential contribution of mineral dust to CDNC, based 
on a simulation that excludes anthropogenic emissions, is up 
to 88.8% over regions of high dust concentrations. Coating 
of dust by hygroscopic aerosol species during its aging 
enhances its CCN activity; a sensitivity study, assuming that 
dust is coated with 10% ammonium sulfate (in volume), 
results in an increase of cloud droplets activated from dust 
particles (by more than twofold). This enhanced hygro- 
scopicity leads to significant water uptake in the early stages 
of cloud formation, resulting in a significant depletion of 
5max (up to 80% over arid areas and up to 30% over their 
surrounding areas). This large depletion also results in a 
significant decrease of total CDNC (up to 70 and 40% over 
the Gobi and Sahara deserts, respectively) but also over their 
surrounding areas such as Southern Europe (up to 5%), the 
Mediterranean Sea, and the tropical Atlantic Ocean (up to 
20%). These numbers represent the first order tendency of 
CDNC changes; assessment of the total response of the 
coupled system would require an interactive aerosol-cloud- 
climate model and is left for a future study. Despite the 
sensitivities and uncertainties seen here, this study shows 
that dust-cloud interactions can be important. Neglecting 
adsorption in predicting the dust CCN activity may intro- 
duce a significant bias in the models that treat dust-cloud 
interactions. Given that dust may affect precipitation in 
climate-sensitive areas, the ability to capture the complex 
impact of mineral dust on cloud droplet formation is an 
important issue for global and regional models. This study 
demonstrates that a comprehensive treatment of the inherent 
hydrophilicity from adsorption and acquired hygroscopicity 
from soluble salts in dust particles is straightforward and well 
within the capabilities of current model frameworks. 
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